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Abstract 

The hybrid spectral problem where the field satisfies Dirichlet condi- 
tions (D) on part of the boundary of the relevant domain and Neumann 
(N) on the remainder is discussed in simple terms. A conjecture for the 
Ci coefficient is presented and the conformal determinant on a 2-disc, 
where the D and regions are semi-circles, is derived. Comments on 
higher coefficients are made. 

A separable second order hemisphere hybrid problem is introduced 
that involves Robin boundary conditions and leads to logarithmic 
terms in the heat-kernel expansion which are evaluated explicitly. 
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1. Introduction 



The explicit construction of the general form of the heat-kernel expansion co- 
efficients has reached the stage when further progress is impeded mainly by un- 
gainliness. Unless there is some compelling reason for finding a specific higher 
coefficient, its exhibition is not particularly enlightening and is not really worth 
the, often considerable, effort. Other, more productive, avenues consist of gener- 
alising the differential operator, the manifold or the boundary conditions. In the 
latter context a simply stated extension is the class of problems where the field 
satisfies Dirichlet conditions (D), say, on part of the boundary and Neumann (N) 
on the remainder. These boundary conditions are sometimes termed 'mixed ' in 
the classical literature {e.g. Sneddon [1]) or sometimes 'hybrid' {e.g. Treves [2] 
chap. 37). A brief history of the corresponding potential theory, sometimes referred 
to as the Zaremba problem, is contained in Azzam and Kreysig, [3]. It is also in- 
teresting to note that these conditions have occurred in string theory, [4] , and have 
recently been considered in connection with isospectrality, Jakobson et al, [5]. 

To set some notation, the conventional short-time expansion of the integrated 
heat-kernel of a smooth boundary value problem is. 



where P is a smooth (singularity free) elliptic differential operator and where, ini- 
tially, the manifold, its boundary and any boundary conditions are all smooth. 
Typically P is the Laplacian, plus possibly a smooth potential, and the coefficients 
are locally computed as integrals over the manifold, or its boundary, of local ge- 
ometrical invariants constructed from the curvature, for example. See [6] for an 
extensive treatment. 

Any relaxation of smoothness can result in modifications to this expansion. For 
example, a singular potential can lead to 'anomalous' powers of t, e.g. [7]. 

The hybrid N/D case has non-smooth boundary conditions and so can be 
classed as a singular boundary value problem, see e.g. [8,9]. Even though it may be 
geometrically smooth, the codimension-2 region, E, where the D and N conditions 
meet, can be regarded as a conical singularity. The general existence of logarithmic 
terms in the expansion of K{t) for such singularities, and other situations, has been 
analysed for a long time {e.g. Cheeger [10], Briining and Seeley [11], Grubb and 
Seeley [12], Grubb [13,14], Gilkey and Grubb [15], Seeley [16]). 
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Seeley, [17], based on [11], has proved the existence of an asymptotic expansion 
that allows for logarithmic terms but then shows for a particular N/D situation that 
such terms do not appear. A classic case is the simple 7r-wedge where the explicit 
calculation shows that logarithmic terms are absent. 

It is anticipated that the heat-kernel coefficients will receive contributions from 
E. This has been confirmed by Avramidi [8,9] and work by van den Berg and Gilkey, 
[18], on heat content is also pertinent. 

In this paper I wish to discuss some aspects of the hybrid question that are 
mainly example driven and with a minimum of algebra. It is hoped that these 
considerations will prove useful in more general field and string theoretic areas where 
heat-kernel coefficients play important roles in divergence and scaling questions. My 
approach is mostly global, as opposed to the local treatment by Avramidi, [8,9]. 

I begin with the Laplacian eigenproblem on the interval with D and N con- 
ditions on the ends. This is then embedded in higher dimensions and used to 
determine the hybrid coeflBcient for the 2-wedge from the purely D, or purely N, 
expressions, which are very old. This result is then used to write down the general 
Ci for a d-manifold with piecewise smooth boundary and a conjecture is made for 
the case where N is replaced by Robin conditions, denoted (R). A crude check in 
the N-D case results from applying the technique to the 2-lune. 

Although my main attention is directed at C^, some very limited information 
on the higher coefficients, C3/2 and C2, is extracted in section 5 from the hybrid 
half-disc. 

As an example of the use of the conjectured form of the hybrid C^, I evaluate 
the Laplacian functional determinant on the N-D disc in section 6. 

In section 7, I set up a separable Robin hybrid problem for the Laplacian on the 
2-sphere and show in later sections that logarithmic terms appear in the expansion 
of the heat-kernel. Perturbation theory is used to bolster confidence in the existence 
of the model. 

2. Basic idea 

For the Laplacian, simple calculation, or the drawing of a few modes, shows 
that on the interval of length L with Dirichlet (D) and Neumann (N) conditions, 
the spectral data of the various eigenproblems are related by, 

{D,N)lIJ{N,N)l^{N,N)^l 
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{D,D)^1}{N,N)l = P^l, (3) 

where the notation (Z), A^) signifies a problem with D conditions at one end and N 
at the other and P stands for periodic conditions. Averaging (2) gives, using (3) 

( A iV) U ^P2L = ^^4L • (4) 

These relations can be extended algebraically to any spectral quantity, such as the 
heat-kernel. 

The 'subtraction' implied by (2) and (4), in order to extract the {D,N) part, 
amounts to a cull of the even modes on the doubled interval, as is well known (c/ 
Rayleigh [19], vol I, p.247). 

The relations can be applied to the arc of a circle, which might form part 
of an SO (2) foliation of a two-dimensional region (or the projection of a higher 
dimensional region onto two dimensions). A wedge is a good example which I will 
now look at. Separability of the modes implies that the relations (2) apply equally 
well to the wedge, where the notation signifies that either D or N holds on the 
straight sides, (say ^ = and 9 = jS). Equation (2) can then immediately be 
applied to the heat-kernel, and its small-time expansion, to determine the form of 
the heat-kernel coefficients in the {D, N) combination. I will show how this works 
out for the coeflBcient. 

The (L>, D) and {N, N) wedge coefficients are well known and have been derived 
in several ways. They are, 

,2 
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^wedge^^^ L)) = Cr^^^N, N) = ^~ ^ , (5) 



Hence from (2) 



cr'^^(Aiv) = -^^^^. (6) 

This last result has been derived by Watson in a rather complicated way using the 
modes directly, [20]. 

Incidentally the conjecture by Gottlieb (equn.(3.5) in [21]), that the {N,D) 
case differs from the (D, D) one only by a sign, is incorrect, although it is true in 
the special case of a right-angled wedge, as is easily checked by looking at rectilinear 
flat domains. The statement is carried through into ref [22]. 

Sommerfeld, [23] vol.2 p. 827, also mentions the 'mixed' wedge and indicates 
how to treat it using images if /3 = nn/m. 
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It is useful to note that, as pointed out by Cheeger, [10] p. 605, the expressions 
(5) and (6), are not locally computable geometric invariants, as evidenced by the 
1/(3 dependence. In the derivation of the wedge coefficients by Cheeger (see Bordag 
et al, [24], Cognola and Zerbini, [25]) the terms arise from the /3-interval ^- 
function evaluated at the argument —1/2. As noted in [24] this is the Casimir 
energy on the interval, clarifying the nonlocal character of this term. By contrast, 
the term proportional to f3 is locally computable. 

It is an important technical point that, as discussed by Avramidi, [8,9], and 
Seeley, [17], it is necessary to specify boundary conditions at the singular region E to 
give a well defined problem. The derivations of the expressions (5) assume Dirichlet 
at the wedge apex. This amounts to taking the Friedrichs extension by default and 
extends to the hybrid wedge (6) by (2). I will continue with this simplifying choice 
for the rest of this paper. 

3. The general case 

Consider in general dimension a manifold whose boundary is piecewise smooth 
consisting of domains, dA4j^, which intersect in codimension-2 manifolds, X^^ . On 
each of the pieces, either D or N is imposed. Kac's principle of not feeling 

the boundary (a cheap way of avoiding estimates), [26], implies that, will take 
contributions from the manifolds of codimension zero, one and two independently. 

In general dimension, for all D or all N, the smeared coefficients are known, 

where dA4 is the union of the boundary pieces and I that of the intersections, 

i i<j 

The 'smeared' coefficients result from the trace Tr /e~^*, where / is a spatially 
local operator, and are handy when discussing conformal transformations and for 
extracting the integrands. I will not use this freedom in any overt calculational way 
(see Kirsten, [27]). 
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For a mixture of D and N, the volume contribution clearly remains unchanged 
while the surface contribution divides simply into a sum separately over those re- 
gions dM.{D) and dM.{N) subject to D and N respectively. The codimension-2 
intersections 2^^- divide into the three (wedge) types X(-D, -D), 1{N^ N) and T{N^ D) 
so the corresponding is, using (6), 

C^{D,N)=(\-di I RfdV+ [ (}^f,-ln-d)fdA 

\0 / Jm JdM{D) \^ ^ / 



JdM{N) V-J ^ / Jx{D,D)UX{N,N) P 

Jj(D,N) P 



In accordance with a previous remark, the wedge-like codimension-2 contribu- 
tions in the above expressions are not locally computable. 

If the boundary is smooth, then all the dihedral angles (3 equal tt and the 
codimension-2 part of (9) (the last two integrals) reduces to 

-\ I fdL. (10) 

4 Jx{D,N) 

To repeat, even though the boundary is smooth, the region T{D, N) = T, is a 
singular region. 

For example, for the 3-ball with D on the northern hemisphere and N {S = 0) 
on the southern, 

for the smearing function, / = 1. 

A local derivation of (10), justifying the application of Kac's principle, has been 
given by Avramidi [8,9]. It has also been obtained by van den Berg (unpublished). 

For later reference, I would like to extend the Neumann conditions to Robin, 

(R), ones. 



{n-d-S)^ 



= 0, (11) 

dM 



where n is the inward normal and S can depend on position. For example, in place 
of (8) one might expect. 
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The first two terms are the standard ones, e.g. [28,6], for a smooth boundary. 
The last, codimension-2, term has actually not been derived directly for the Robin 
wedge but it holds when S = and dimensions show that the Robin function cannot 
enter algebraically into this term. 

On the same basis my conjecture for the corresponding C^{D, R) is, 

C,{D,R)=(l-^) f RfdV+f (l^-L-d)fdA 
J Jm JdM{D) \^ ^ J 

+ [ (^^-2S + ln-d]fdA + l [ t^jdL 

JdM{R) \^ ^ J ^ JX{D,D)\JI{R,R) P 

7r2 + 2,?2 



12 Jx(D 



-fdL. 



:(D,R) P 

(13) 

Proceeding on the basis that (13) is correct, an expression for Ci{N, R) can be 
obtained by making dM.{D) empty and dividing dM.{R) into dM.{N) U dM.{R). 
the conjecture is then, 

^ (14) 



C,(N,R)=(l-()I^Rf4V 

JdM(R) ^ Jx P 



'dM{R) 

It is clear that the above expressions have a specific validity, even without the 
codimension-2 parts. Thus, although the Robin form (12) trivially reduces to the 
Neumann one, (8), when S' = 0, it is not possible to obtain the Dirichlet form, (7), 
directly by setting S = 00. The exhibited forms, which, to repeat, are coefficients 
in a 'small time' asymptotic expansion, are really valid in the limit of small S'^t, 
as discussed by Fulling, [29]. I return more specifically to Robin conditions from 
section 7 onwards 



4. The lune 



The expression for Ci{D, N) can be checked in a curved space case by consid- 
ering a lune segment of a sphere. 

Relations (2), (3) can be applied to the 2-lune where the intervals are the 
sections of the lines of latitude cut out by the the two longitudes, (/> = 0, (/) = (3. In 
this case the extrinsic curvatures vanish (the boundaries are geodesically embedded) 
but there is a volume (area) term independent of the boundary conditions. 
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The ^-functions are now somewhat more exphcit [30,31]. It is possible to work 
with general angle (3. If one chooses /? = tt/q', e Z, the ^-functions have been 
derived in [31] and used in [32]. 

Denoting the lune by C{(3) one has, 

(AiV)z;(;3) U (A^)r(;3) = (AI^)z:(2/3) 
(A iV)z:(;3) U (iV, iV)^(^) = (iV, iV)^(2^) , 

so that the corresponding (^-functions combine algebraically, 

cr(-) - C2T(-) - cr(-) - C2T(-) - cr(-) • (15) 

The DD and A^'A'" ^-functions have been derived in [31] as Barnes ^-functions 
for conformal coupling in three dimensions (leading to simple eigenvalues) and yield 
the specific value, for example. 



which can be used to confirm the expression (6) using the relation between C-^ and 
C(0). (In this case there are no zero modes.) 

The volume contribution, /3/6, to is standard and is the same for all bound- 
ary conditions. Hence the contribution of each [N^ D) corner (of which there are 
two) is 



1 
2 

as required for the check. 



47r/7r /3\ /? 



24V/3^7ry/ 6 



12/3 



5. The disc and semi-circle. Higher coefficients. 

The fact that the extrinsic curvatures arc zero means that the lune is not 
excessively helpful in deriving the form of the higher coefficients in the (D, A^) case. 
Some further information can be obtained by looking at the half-disc with semi- 
circular boundary having different conditions on the diameter and circumference. 

A straightforward application of, say the Stewartson and Waechter Laplace 
transform technique combined with an image method soon yields the results for the 
short time expansions 

, ^ 1 2-f-7r 5 v^(7r + 16) f I 1\ 
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r ^ 1 2-7r 1 v^(7r-16) /I l\ 

/ N 1 2-7r 1 v^(57r-48) /I 3\ 

where DA?^ means on the diameter and N on the circumference, etc. 

The constant terms check with (5) and (6) for (3 = 7r/2. Also (3), applied to the 
diameter as a wedge of angle tt, yields the D and N (e.g. [33]), full disc expansions. 

The extrinsic curvature vanishes on the diameter and equals one on the circum- 
ference part of the boundary so some information on the C3/2 and C2 coefficients 
can be inferred. Formulae in the non-mixed types {D,D) and {N,N) have been 
given in [32,34] which agree with the relevant parts of the above expressions. Indeed 
I used the hemi-disc in deriving these results. 

Also in [32] will be found an expression for C2 in the case the boundary parts 
dA4j^ are subject to Robin conditions with different boundary functions, 5^ although 
all dihedral angles are restricted to 7r/2. 

In the case of C2, the 1/315 is the contribution of the curved D semicircle while 
the ±1/32 is the effect of the (two) corners and likewise regarding the 1/45 ± 3/32 
combination. The C3/2 coefficient exhibits a similar structure. Experience with the 
flat wedge shows that it is unwise to draw too many conclusions when the angle is 
7r/2. What we can say, however, is that, using the 3/2 coeflicient as an exemplar, 
one term will have the general form 

/ ^DD{l3){i^i + i^2)dL+ j Xj^^{P)(k-^ + K2)dL+ 

Jx{D,D) Jx(N,N) 

{^nd{(^)i^d + >^dn{(^)i^n) dL 



24 



I 



T{N,D) 

where A^o(7r/2) - -A^^(7r/2) = 3 and A^^(7r/2) = -A^^(7r/2) = -9. This 
change of sign is a simple consequence of images, or of (2) since the DD and A^'A'" 
quantities vanish when j3 = tt. 
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6. The disc determinant. 



A direct attack via modes, of what is, after rectilinear domains, the simplest 
two-dimensional situation, i.e. a disc subject to on one half of the circumference, 
and D on the rest, would seem to be difficult in so far as the construction of the 
function or heat-kernel is concerned. However, the functional determinant, defined 
conventionally as exp ( — C'(0)); appears to be accessible by conformal transforma- 
tion from that on an ND-\\ine of angle tt, i.e. a hemisphere with N on one half of 
the rim (the equator) and D on the rest, which is an easy quantity to find in terms 
of Barnes ^-f unction from (15). 

Instead of the determinant I use the effective action, W, defined hy W = 
— C'(0)/2. Integrating the conformal anomaly leads to the relation 

W[e-^^g]=W[g] + W[e-^^g,g] 

where W[e~'^'^ g,g] is the cocyle function. 

For the above programme to work, one would need the conjectured form of Cj^, 
(13), to be valid in order to construct the required cocycle function in two dimen- 
sions. Applying the standard techniques this is (c/ [35]), for a smooth boundary, 

W[e-^^g,g]^^ f u;{R + nu) dV + [ iv{K + l{n.d)u;) dA+ 

2471 Jm 127r Jqm 2 ^^^^ 

^( f -f ]{n.d)ujdA-^y2uk 

°'^\JdM(N) JdM(D)/ u 



where k labels the points where D and A^ meet and are the values of uj at these 
points. If dAi{D) is empty there is a volume term coming from the pure A^ zero 
mode. 

To go from the hemisphere to the disc I employ the equatorial stereographic 
projection as in [36,37,32,34,35] noting that there is no codimension-2 contribution 
because the conformal factor is unity on the boundary, implying u;^ = 0. 

Then (20) can be written 

W:,n[9,9] = l(wn[9,9] +Wr,[g,g]), (21) 



where VFjv means the usual Neumann expression, omitting the zero mode piece, 
and I can use the known values, (g = disc and g = hemisphere), 

WD[g,9] = l^og2-\ 
^Ar[^,^] = ^log2+^. 
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The ^-function on the A'"D-hemisphere follows from (15) with (3 = tt. The 
(■-function, C,^^{s) is the usual hemisphere ^-function and the determinant has 
been considered a number of times. C,2tF{s)^ corresponds to Sommerfeld's double 
covering of three-space introduced in connection with the half- plane boundary. 

Since one needs conformal invariance in two dimensions, not three, the (- 
functions are actually modified Barnes (-functions which have been dealt with in 
[30,38,39]. The determinants can be computed generally in terms of Barnes (- 
functions but, because of the rational nature of 7r//3, in this case, they can be 
reduced to Epstein or Hurwitz (-functions. The general theory, appropriate to 
the arbitrary 2-lune, is developed in [30]. However it is probably easier to proceed 
directly. 

From [30] the (-function for —A on the ND 2-hemisphere is 

oo ^ 

/■ND( \ _/-DDf \ /-D-D/ \ _ \^ 

((1 + m + n)2 - 1/4) 
Expanding in the 1/4 leads to the expression for the derivative at 0, 

C'(0) = C^(0, 1/2 I 1, 1) + (^(0, 3/2 I 1, 1) - ^ , (23) 



where 



oo ^ 

Cn(s, a \ 1,1)= > 7 r- 



m,n=0 

is a 2-dimensional Barnes (-function and N2{a) is its residue at s = 2; N2{a) = 1. 
In this simple case the sums can easily be rearranged, 
oo ^ °° N +1 

^ (a + m + nY ^ (N + aY ^ shv , y , \ j 

m,n=0 ^ ' N=0 ^ ' 

SO that 

(2(5, 1/2 I 1, 1) + (2(5, 3/2 I 1, 1) = 2(^(s - 1, 1/2) = 2(2^-1 - 1) (^(s - 1) 
and therefore from (23) 

(f^'(0) = -(^(-l)-^log2-i. (25) 

The absence of a (/^(O) term is related to the absence of the perimeter heat-kernel 
coefficient caused by the equal-sized and D regions. The (-function has only the 
Weyl volume pole. 
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For comparison the standard formulae for the DD and A^'AT-hemispheres are 

cr'(0) = 2C^(-i)-CM-^ 

and 

Cf^'(o) = 2C^(-i) + ao)-^. 

By conformal transformation, on the A'"D-disc, the final result is 

rdisc 1 , -X 11 , , 1 



<lf = 5a-i) + 5jiog2-^, 



using (21) with (22). 



7. The Robin boundary condition. f 

The Robin condition, (11), has made only a formal appearance in the discussion 
so far. It was needed for conformal transformations but has not yet entered into 
any eigenproblem. 

The reason why Robin conditions are so awkward is that, in general, the prob- 
lem does not separate and, even if it does, the eigenvalues are given only intrinsically. 
Early considerations of the eigenproblem are reviewed by Pockels, [43]. Poincare, 
[44], also used the condition in connection with eigenfunction existence. See also 
Bandle, [45]. A practical, more recent, treatment is given by Strauss, [46]. 

Apart from applied mathematics, there has been some recent interest in Robin 
conditions in the quantum field theoretic and spectral geometry scenes { e.g. Fulling, 
[29], Bondurant and Fulling [47], Romeo and Saharian [48], Solodukhin, [49], de 
Albuquerque and Cavalcanti, [50]). 

In this section I return to the 2-hemisphere, some aspects of which were men- 
tioned in the previous section. Here I wish to see how far the standard Robin 
eigenproblem on the interval, {e.g. [43,46,51]), is relevant for a spherical geometry. 

Appropriate details of the classic hemisphere eigenproblem were also given in 
[52]. To repeat, cooordinates on the hemisphere are < ^ < tt and < < tt. 
The rim (boundary) consists of the union of the two semicircles = and = tt. 
In [52], D conditions were applied at = 0, and N at </> = tt. Now the latter are 



t According to Gustafson, [40,41], Gustave Robin (1855-1897), never seemed to have used this 
condition. His name was, apparently, first attached to it by Bergmann and Schiffer in the 1950's 
but the condition occurred already in the work of Newton [42] . 
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replaced by Robin conditions and the former by either D or N. It is also possible 
to treat both Robin, but I will not make this generalisation simply for convenience. 
The singular region, E, comprises just the S and N poles. 
The Robin condition (11) specifically is 



1 a$ 



(26) 



sin 6 dcj) 

which is not consistent with a separated structure for ^> unless S takes the form 

^ = -^ (27) 
sm ^ 

diverging on E. I will return to this point later and take S as in (27) simply in 
order to get on with the calculation because, in this case, the condition (26) reduces 
to the usual (D,R) (or (N,R)) on the (j) 'interval' (0, tt) and I can employ known 
results. In the separated solution for $ the 9 part is unchanged, only the 4> factor 
is modified. Thus, in the (D, R) case, the hemisphere eigensolution is 

$;,=A/;sin(A:0)p-+\(cos^), (28) 

where A; > is determined by the one-dimension interval condition, 

k(X)t{k7r) = h, (29) 

so that, as /i — 0, /c tends to half an odd integer, k —>■ m + 1/2, m = 0, 1,..., the 
Neumann result. 

Likewise, for the {N, R) case, instead of (28) and (29) there is 

^^=X^cosik<j>)p-^,{cose) (30) 

and 

ktan{k7r) = -h. (31) 

This time, as /i — > 0, /c — > m, m = 0, 1, . . .. In both cases I will label k by the 
associated m. 

To be specific, the eigenproblem I now wish to consider is 

(_A + i)$, = A$;,, (32) 



12 



the reason being that the eigenvalues are perfect squares, 

1 2 

^ = ^mn= + , m,n= 0,1,.... (33) 

Degeneracies, if there are any, are due to coincidences. 

For (A^, R) conditions with /i > 0, /cq is imaginary but does not correspond to a 
hemisphere mode because the assumed self-adjointness imphes that the eigenvalues 
A must be real. Then, in this particular case, m starts at 1. 

One aim is to relate the spectral properties of the A to those of the by 
summing out the n. I will do this via the heat and cylinder kernels. We have used 
this before in spherical situations, [31]. The heat-kernel and cylinder (or 'square 
root') kernel are defined, in general, by 



xt 



(34) 

A 

where I have included a degeneracy, just in case. From now on I use FuUing's 
notation, setting T{t) = i^^/^(t), and taking t as a generic parameter. 

Using the expression (33), it readily turns out that the hemisphere cylinder 
kernels factorise, 

Trrc;(t) = ^——Tr{t) , (35) 

^ 2sinht/2 ^ ' ^ ' 

Tj being the cylinder kernel on the interval defined by, 

oo 

Ti{t)= (36) 

m=mo 

where uiq = 1 for (AT, R) with h > and uiq = otherwise. Equation (35) is the 
connection between the hemisphere and the 7r-interval. 

Important information is contained in the short-time expansions of the heat 

1/2 

and cylinder kernels. A reflection of the pseudo-operator character of ( — A +1/4) 
is the possible existence of logarithmic terms in the expansion of Tjjg. 

It is relatively straightforward to show that on a d-dimensional manifold, 
maybe with a boundary, the expansion of a general cylinder kernel, T, takes the 
form, 

oo oo 

T{t) ~ J] a, + t-'^' log t (37) 

i=0 i 
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The lower limit on the second term is deliberately left unspecified but I draw 
attention to the important fact that, if the operator P in (34) is a smooth differential 
operator with smooth boundary conditions, then only odd positive powers of t occur 
in the log term, whatever the dimension of the manifold. 

One way of showing this is to use the known existence, in this case, of the series 
expansion of the heat-kernel K{t), (1), and then, by (^-function manipulations, 
relate the coefficients a, a' and C. This was first done in a physical context by 
Cognola et al [53] . A more recent analysis, in the compact case, has been performed 
by Bar and Moroianu [54] who consider local, diagonal kernels and give a careful 
analysis of estimates. Another way is to relate the asymptotic expansions obtained 
by smoothing using either A or v^A as the preferred variable. This was employed by 
Fulling [29]. It is sufficient to note, for now, that the coefficients of the logarithmic 
terms, a[ in (37) are determined by the hj^. 

The heat-kernel coefficients, 6^, for Kj, on the {R,R) interval have been ob- 
tained in [55] and show that there are logarithmic terms in Tj. The relation (35) 
between cylinder kernels then implies that Tjjg also has logarithmic terms, but with 
even powers of t. In turn, this suggests that the operator —A + 1/4, together with 
the boundary conditions, including the choice of Robin parameter, (27), is a rather 
singular operator. This is looked at further in the next section. 

8. Asymptotic series. 

To allow for the fact that the operator — A + 1/4 might be particularly singular, 
I generalise (1) to include logarithmic terms, the immediate aim being to relate the 
expansions of the heat and cylinder kernels. For this purpose, I use the zeta-function 
approach, mentioned earlier, together with the general series established by Grubb 
and Seeley [12]. Suitable summaries can be found in [56,57]. For notational brevity 
I set P = —A + 1/4 and Q = \fP ^ but, so far as the general equations go, P can 
be any Laplace-like (second order elliptic) operator of smooth form. 

The asymptotic expansion of the heat-kernel i^(t), (34), is, e.g. [56,57], on a 
d-dimensional manifold, 

oo 

-d<fc<-d+2 k=-d-V1 

The reason for the lower limit of —d + 2 will appear later. This limit differs from 
that in [11,57], which is zero. 
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The simple powers have been spht into two sets because the coefficients have 
different qualities. Since this does not concern me at this time, I will combine them 
for algebraic ease. Therefore, 

oo oo 

m^j2ht'/'+ E b^f^^'^^st. (38) 

k=-d k=-d+2 

This generalises (1), with the relation between the coefficients 

Similarly for the cylinder kernel {Q is a first order pseudo-operator) I will assume, 

oo oo 

T(t)~ J2''kt''+ E 4*'logt, (40) 

k=-d k=-d+2 

where now, all powers of t greater than —d + 1 appear in the logarithmic term. 
The connection is made via the corresponding zeta-functions, 

Cq(2.) = Cp(s), (41) 

which have asymptotic expansions corresponding to (38) and (40). I refer to [56], 
e.g. , in order to save work, and deduce, 

oo , oo , ^ 

oo oo ; 

k=-d k=-d+2 ^ ' 

where ng is the number of zero modes. 

The relation between the a's and the 6's follows from (41) using the standard 
duplication formula for r(2s). Replacing s by 2s in (43) it is easy to see that dividing 
by r(s + l/2) removes certain first order poles and converts some second order poles 
into first order ones. From the residues, making the necessary identifications, I find 
the relations 

2^ I — 

^k = —r, TT^ (^ki ^ = • • • — 1, 0, 2, 4, . . . 

^ r((l-A;)/2) ^ 



2"-^v/^ 

r((i-fc)/2) 

/c = l,3, 



(44) 



^fc = TTT:^^' /c= -d+2,...,-l,0,2. 
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which generahse those found by FuUing [29] and Cognola et al, [53] . 

Some overall conclusions can be drawn from these relations. Firstly, given the 
heat-kernel bj^, it is not possible to determine the cylinder a^. for odd positive, 
k = 1,3,..., as emphasised by Fulling. Secondly, for the assumed structure of the 
cylinder expansion (40), or (43), one cannot specify all the logarithmic terms in 
the heat-kernel, i.e. all the 6'^. Since the form (40) is sufficient for the quantities 
appearing in this paper, I will leave this point except to say that it is easy to take 
(42) and work equation (41) the other way to derive the corresponding expansion 
for Cqi^)- The new feature is the appearance of poles of third order, leading to 
higher logarithmic powers in the cylinder kernel expansion. 

The extension of the lower limit down to —d + 2 has the rather violent conse- 
quence that, in this case, the coefficients 6j for i > —d + 2 are global {i.e. nonlocal) 
quantities. Only b_^ and are locally computable. These are the Weyl volume 

and boundary area terms denoted by Cq and earlier. This means, in particular 
that Ci, discussed in §3, is not locally computable. We have seen that this is true 
for the exhibited form, (13), through the last two, i.e. codimension 2, terms. The 
question of nonlocal terms is considered in a later section. 

A further consequence of (44), and the one relevant for the spherical problem 
treated in the previous section, is that the logarithmic terms in the heat-kernel arise 
from those logarithmic terms in the cylinder kernel with even powers of t, which is 
precisely the case with (35). 

The conclusion is that the Laplace operator, —A + 1/4 with the boundary 
conditions, is a rather singular operator as it provides a concrete example of a 
second order problem involving logarithmic terms. It would therefore seem that it 
does not come within the compass of Seeley's analysis, [17], probably because of 
the divergence in the Robin function, S. Hence, before giving explicit forms for 
the expansions, I break off to look at the modes, (28), a little more closely, one 
reason being that, although self-adjointness depends on the formal subtraction of 
two (identical) terms, 

de^S-^, (45) 

at the boundary, the divergence of S, (27), at the poles {6 = 0, tt), i.e. on E, might 
give one pause for thought. In general terms, the assumed Dirichlet conditions at 
E are actually sufficient for convergence of the integral, (45). 
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9. Robin mode properties. 



The orthonormality of (28), and of (30), is easily established either by direct 
integration or, formally, by the usual self-adjoint Liouville method. 

For convenience I write down some standard things. For the Legendre functions 
orthonormality reads 

C dxP-^ ix\p-^ (x')- '2''^in + l) 

ax i^^+^Kx) - ^ 2„ ^ ^^^^2^ + n + 1) ■ 



nn' 1 



where n and n' are positive integers or zero and k > —1 (MacRobert [58] p. 335). 
I also note the limiting behaviour at the poles, 

The interval Robin modes are standard, {e.g. Carslaw and Jaeger, [51], Strauss 
[46], Pockels, [43]) with easily determined normalisations. I again remark that 
Dirichlet conditions continue to apply on E, my default position. 

Use of the limit (46) shows that each boundary term, (45), in the self-adjoint 
condition applied to two eigenfunctions, corresponding to k and /c', is finite if k+k' > 
0. So we are completely safe in this case. 

This limit also implies that the total heat flux, per mode, 

dM 

is finite. Moreover, the 'quantum mechanical flux', 

'dM 

is also flnite, by virtue of Barnes' formula, 

',,M!_±ni±ii±4, (47) 

1-x^ 2// r(l - + z/) ' ^ ^ 

valid for Re n < and fi + u a positive integer, or zero. The lower limit can be 
extended to —1 using the fact that the Legendre functions are unchanged, up to a 
sign, under x —x ([58], p. 334 Examples (2) and (3)). 

Confidence in the eigenmodes is increased if one employs a perturbation tech- 
nique to calculate the change in the eigenvalues of (33) for a small change in 
S. The formula, which will not be derived here, is 



5X= / ^A-^^A- (48) 

JdM{N) 
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For simplicity, I am considering the (D, R) set-up and am perturbing about the 
{D, N) case, so S is small. The integration over the boundary encompasses only the 
part on which R (equivalently (AT)) holds, since ^» is zero on the D part. Equation 
(48) seems to occur first, for constant S, in Poincare [44] and for variable S in 
Pockels, [43] p. 178, as a quick consequence of Green's formula. See also Frohlich, 
[59], equn. (6d). 

The answer is known from a direct analysis of the interval eigenvalue condition 
(29) which shows that 

1 2/i 

2 [Am + l)7r 

Alternatively, applying (48) yields 



(50) 



where m = 2m + 1. The normalisation is, 

2 _ 2 r(m + n+ l)r(m + 2n + 1) 
~ n 2r(n + l) ' 

and substitution of (47) into (50) easily produces 

A /Y _ 2/1 1 

which is just (49). This limited check inspires a certain confidence in the sensibility 
of the modes and the model. 

I should point out now that I do not attach any serious significance to the 
model. The choice S is simply one of convenience for solvability. 



10. Explicit expressions 

In this section I use the form of the heat-kernel coefficients for the interval 
(-R, R) Robin problem derived in [55], from which one can easily deduce those for 
the (-D, R) and {N, R) cases by appropriate limits. 

I find, for both {N, R) {h < 0) and (L>, R) on an interval of length tt, (I), the 
heat-kernel coefficients. 
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so that, for example, the coefficients of the logarithmic terms in the interval cylinder 
kernel, (36), are determined from (44) applied to the interval, as 

and these odd terms are all there are because there are no logarithmic terms in the 
interval heat-kernel. The asymptotic series can be summed to produce the closed 
form for the asymptotic logarithmic part of Tj, 

e-'^'*' - 1 

and then from (35) the logarithmic part of the {N,R) {h < 0), or [D,R), hemi- 
sphere, (HS), cylinder kernel follows as, 

Th's it) ~ ■ w — - log t . (52) 
^ sinh(t/2) 27Tht ^ ^ ^ 

The same result also holds for the {N, R), h > 0, case even after dropping the 
imaginary /cg interval-mode. The contribution of this to the heat-kernel coefficients 
must be removed but this does not affect the 6^, for odd k. 

The relation (44) can now be applied to the hemisphere and the logarithmic 
terms in the heat-kernel worked out. I do not give the easily derived expressions. 
They are combinations of Bernoulli numbers. Their existence is all that is required 
for now. I only remark that the important term ~ t°logt is present. 

There are, of course, a number of technical routes to the expressions and con- 
clusions derived above. I have chosen to use heat and cylinder kernels but one 
could employ the resolvent and a standard contour method of rewriting eigenvalue 
sums. I also note that it is straightforward to extend the calculations to a lune and 
also to d-dimensions, which just gives higher powers of sinh(t/2) in (52), say. This 
justifies the lower limit in (40) and hence in (38). Odd spheres are not singular, in 
agreement with a general result. The details will be presented at another time. 

11. Nonlocal terms and the Casimir energy 

In this section I enlarge on previous statements regarding the nonlocality of 
some expansion coefficients, in particular of which, in two dimensions, is the 
coefficient of the constant term in the heat-kernel expansion. 
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As a first step, I look at the Casimir energy, E, on the 7r-interval (I) and simply 
quote the formal definition. 



E = FP^Cp(-1/2), (53) 

where P = —d?/d(fP' and the boundary conditions are either (D,R) or (N,R) on the 
ends. FP stands for 'finite part'. The Casimir energy is a nonlocal quantity. 
On the TT-interval the relevant definition here of the ^-function is, 

m=mQ 

For the purposes of this paper, I consider E, defined by (53), simply as a conve- 
nient mathematical quantity rather than as something having physical, operational 
significance. In particular, on the pure interval, there would be the question of what 
to do with the 'negative mode' that occurs for (A^, R) when h > 0. In the present 
model, this has already been excluded. (Compare with reference [49] where these 
modes are referred to as 'bound states'.) 

Equations (41), (42) and (43) hold for the interval. Since P is smooth, there 
are no log terms in the heat-kernel and so the coefficients 6/ are zero. 

It is important to realise that in the Robin case, Cp{s) has a pole at s = — 1/2, 

(^iJs) ^ + B (55) 

where the residue, A = — 6{/2y^, follows from (42) and the remainder, B, equals 
2E, by definition. In terms of the interval heat-kernel coefficients, C^, (see (51), 
(39)), 

for both {D, R) and {N, R), for all h. 

From (53) and (41) it is required to work around s ~ — 1 for Cq(^) where 

and so 



+ 7-1 -77^ 777T + ^ 



s + 1 ' J\s/2 + l/2 J s + 1 (s + 1) 
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which yields 

6f = -V^a'/, (58) 

and also 

Equation (58) agrees with (44) while (59) relates the coefficient a{ to the nonlocal 
Casimir energy. It differs from the result in [29] by the last (constant) term. See 
[53] for a relevant discussion of various regularisation recipes in this context. 

The next step is to relate the interval and hemisphere expansions using (35). 
Expansion of the 1 / sinh easily gives the connection, 

«r = 4- , (60) 

which I now rewrite in terms of heat-kernel expansion coefficients. 
Application of (44) to the two-hemisphere gives, 

Qq — Oq , 

where b^^ is the constant term coefficient. 

Equation (44), applied to the interval, gives, 

=—b^ =—-^ = 1 

and so, finally, (60) becomes, using (59), 

Cf^ = -SttE- I -4/i(7- 1), (61) 

which relates the coefficient on the hemisphere to the Casimir energy on the 
interval for my model, i.e. (27). 

A basic check sets h = when (61) allows a computation of the interval Casimir 
energies for the (D,D) and (D,N) cases using the expressions for Ci, (7), (12) and 
(13). Simple algebra yields the standard values, 

EiD,D) = EiN,N) = -^, EiD,N) = j^. 
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It is worth noting that for the (D,D) and (N,N) cases, although the Casimir 
energy, E, is nonfocal on the interval, the coefficients are local on the hemisphere. 

In the general case, is a nontrivial function of h on the interval and the 
conjectured forms of Ci{D, R) and Ci{N, R), (13), (8), obviously do not agree with 
(61). However, it will be recalled that the heat-kernel expansion is really one in 
h'^t and therefore one should treat h as 'small'. As mentioned, it is in this realm 
that (13) and (12) should be valid. Furthermore the divergence of the hemisphere 
Robin function, S of (27), makes Ci formally infinite, and cannot be considered 
'small'. Hence it is clear that (13), (12) and (8) are not even appropriate for the 
present situation. Against this must be set the fact that, as described in section 9, 
perturbation in S appears to work for the eigenvalues. Therefore I intend to give a 
further look at perturbation theory. This will give us a simple, if limited, analytical 
handle on the ('-function on the Robin interval which might also be useful in other 
circumstances. 



12. Perturbation approach 

I first consider the (N,R) case when one must distinguish between positive and 
negative Robin parameters, h. If h > one has the approximation for the interval 
frequencies, 

«^m~™ ' m=l,2, 

rmr 

which leads to 

Oho 

Ci>(s)«^Cfl(2s) + — Cfl(2s + 2), (62) 

TT 

in terms of the Riemann C-function, and so A — —h/27v with 

E(iV,i?)«.-l- A(7-l), hiO. (63) 

By contrast, if /i < 0, there is a root which tends to zero as /i — > and becomes 
the zero (N,N) mode. Approximation of (31) gives. 



h 

^m~™ ' m = l,2, 

rmr 



whence A= —h/2'K and 
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Turning to the (D,R) case, to order h, (49), 



m 2h w? 2h 

~ IT - ^ >^^—. > m = l,3,..., 

so that 

9h<t 

CUs) ^ (2^^ - 1) CrC^s) + — (2^«+2 - 1)C^(2. + 2) . (65) 

TT 

Working around s ~ —1/2 yields, in accordance with (55), the residue check, A = 
—h/2TT, and the Casimir energy, 

£;(i^,i?) ^ I ^ - A(7 - 1 + 21og2) . (66) 

From these expressions, I use (61) to compute the heat-kernel coefficients on 
the hemisphere. I find, in the two cases, the values, 

C^{N,R)^^, h>0 

Ci(7V,i?) f« ^ -4V'^, h<0 (67) 
C^{D,R) ^ -| + 8/ilog2. 

It appears that the effect of the Robin function (27) for /i > has disappeared from 
Ci {N, R) leaving just the {N, N) value, at least in this lowest approximation. 
This asymmetry can be traced to the omission of the negative mode ior h > 0. 
(Compare with the discussion in [49].) 

The log 2 term in Ci{D, R) reinforces the conclusion that (13) is inappropriate 
for the present singular model. 

Similar results can be shown to hold for the simple wedge of section 2 with 
Robin conditions on one side. In order for the techniques of [24] to work, separability 
demands that the Robin function is again singular on E (the apex of the wedge) 
taking the form S = —h/r. Without going into details, just referring to equations 
(4.4)-(4.6) of [24], there is again a log term coming from the pole at s = —1/2 in 
the interval Robin (^-function and there is a log 2 term in the expression for the 
corresponding C^. 
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13. Exact form of interval Casimir energy. 



The exact expressions for the Casimir energy derived in [48] and [50] can be 
approximated for comparison with my perturbation results. For convenience I refer 
to equations (34) and (39) in [50] for the Casimir energies and rewrite them in one 
dimension, 

S{N, R) = E{N, N) + — [ dk log (1 - 



£{D, R) = E{D, N) + ^ I dk log 

2tt J 



27: J "V (A;-/i)(exp(2A;7r-l))/ 

2k \ 1 ^^^^ 

1 + 



(A;-/i)(exp(2A;7r-l))y 16' 

where the S may or may not, agree, up to renormalisation, with the quantities 
evaluated from (53). (I have set h = — C2 < 0, in the notation of [50], to comply 
with my sign.) 

The leading small h behaviours of the 'correction' terms in (68) can be deter- 
mined numerically to be. 



S{N,R)^E{N,N) + lJ-^ + ^F^ih) 

2 V TT 27r (69) 

£{D,R)^E{D,N) + ^Fi,{h), 
where Fpj- and F^, satisfy the functional relation 

F{Xx) - XF{x) ^ ^ . (70) 
A 

There is some similarity between (69) and the perturbative results, but I cannot 
explain (70) from such a viewpoint and simply present the result as a possible 
significant curiosity. 



14. Conclusion 



Apart from rectilinear domains, and the hemisphere, there seem few situations 
that can be solved exactly for A^D-conditions (see [43]) and this is a drawback to 
the construction of the precise forms of the heat-kernel coefficients. Nevertheless 
a certain progress has been made in a simple minded way making use of the ND- 
wedge expression. This type of reasoning can be extended to higher dimensions 
leading to information about the trihedral corner contributions and their higher 
analogues. 
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Surprisingly the conformal functional determinant is available on the 'half-N 
half-D' disc by conformal transformation from the A/"Z)-hemisphere, and has been 
computed, assuming that a conjecture for the heat-kernel coefficient, C^, is correct. 

In sections 7 to 10, I considered hybrid {D, R) and {N, R) 2-hemisphere prob- 
lems, with a Laplace style operator, which have logarithmic terms in the short-time 
expansions of the heat-kernels. The boundary conditions involve a Robin function 
that diverges, but not too strongly, at the poles which are the places where the 
D, N conditions meet the R condition. This allows the problem to be separated. 

Although the second order differential operator (with the boundary conditions) 
is rather singular, it is not singular enough to prevent the heat-kernel expansion 
from existing. 

The (singular) Robin condition is essential here for the existence of the loga- 
rithmic terms. The expression (52) vanishes when h = 0. In [52] it was suggested, 
by an indirect argument, that for an {N, D) situation, a nonzero extrinsic curvature 
at the boundary gave rise to logarithmic terms. In general, it seems that the Robin 
condition mimics an extrinsic curvature, [29] equn, (4.1). If this is so, it is perhaps 
not surprising that a singular Robin function, S, produces a log term as it would 
simulate a conical type singularity. 

The attempt to relate the hemisphere heat-kernel coefficient to the Robin 
Casimir energy on the interval, analysed in sections 11 and 12, confirms the lim- 
ited validity of the expression for the coefficient, (13), to extend which requires 
further analysis. 

An incidental line of technical enquiry would centre on the significance of the 
'imaginary modes' on the {N, R) interval for h > and whether continuity could 
be maintained through h = 0. 

Further aspects of this particular set-up will be examined at another time. 
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